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Abstract —Free space optical (FSO) communication has been 
receiving increasing attention in recent years with its ability to 
achieve ultra-high data rates over unlicensed optical spectrum. A 
major performance limiting factor in FSO systems is atmospheric 
turbulence which severely degrades the system performance. To 
address this issue, multiple transmit and/or receive apertures 
can be employed, and the performance can be improved via 
diversity gain. In this paper, we investigate the bit error rate 
(BER) performance of FSO systems with transmit diversity or 
receive diversity with equal gain combining (EGC) over atmo¬ 
spheric turbulence channels described by the Double Generalized 
Gamma (Double GG) distribution. The Double GG distribution, 
recently proposed, generalizes many existing turbulence models 
in a closed-form expression and covers all turbulence conditions. 
Since the distribution function of a sum of Double GG random 
variables (RVs) appears in BER expression, we first derive a 
closed-form upper bound for the distribution of the sum of 
Double GG distributed RVs. A novel union upper bound for 
the average BER as well as corresponding asymptotic expression 
is then derived and evaluated in terms of Meijers G-functions. 

Index Terms —Free-space optical systems, atmospheric tur¬ 
bulence, Double GG distribution, bit error rate, equal gain 
combining, spatial diversity. 

I. Introduction 

F REE-SPACE optical (FSO) communication, which has 
been receiving growing attention, refers to terrestrial line- 
of-sight optical transmission through the atmosphere using 
lasers or light emitting diodes (LEDs). FSO systems offer cost- 
effective, license-free and high capacity communication which 
is a promising solution for the ’’last mile” problem. The unique 
properties of these systems also make them appealing for a 
number of other applications including fiber backup, wireless 
metropolitan area network extensions, broadband access to 
remote or underserved areas, and disaster recovery. EMU- 
A major performance impairing factor in FSO systems is 
the effect of atmospheric turbulence-induced fading resulting 
in random fluctuations in the received signal (4j. Atmospheric 
turbulence occurs as a result of the variations in the refrac¬ 
tive index due to inhomogenities in the temperature and the 
pressure of the atmosphere. In the literature, many statistical 
models have been proposed to model this random phenomenon 


for different degrees of turbulence severity. Under weak tur¬ 
bulence conditions, log-normal distribution has been widely 
used to model the random irradiance fluctuations EMU. As 
the strength of turbulence increases, log-normal distribution 
shows large deviations from the experimental data. Therefore, 
other statistical models such as K ED, i-K EH, log-normal 
Rician lll3ll . Gamma-Gamma fl4l . M 031 and Double Weibull 
OH distributions have been proposed to cover a wide range 
of turbulence conditions. Recently, a unifying statistical distri¬ 
bution named Double Generalized Gamma (Double GG) was 
proposed in fl7l which is valid under all range of turbulence 
conditions and contains most of the existing models in the 
literature as special cases. 

Spatial diversity techniques provide a promising approach 
to mitigate turbulence-induced fading. Over the years, the 
performance of FSO systems with spatial diversity have 
been extensively studied over the most commonly utilized 
turbulence-induced fading models including log-normal, K, 
negative exponential and Gamma-Gamma Ifl8l-ll22l. In 03, 
a unified closed-form expression for the BER of single-input 
multiple-output (SIMO) FSO systems with optimal combining 
(OC) receiver over Double GG channel was proposed. How¬ 
ever, the BER performance of equal gain combining (EGC) 
receivers, which is equivalent to that of multiple-input single¬ 
output (MISO) FSO systems, was presented in integral form. 
In fact, the main difficulty in studying EGC receivers is the 
distribution of the sum of Double GG R.Vs required to be 
derived in order to obtain a closed-form solution. 

An analytical solution for the distribution of the sum of 
Double GG RVs is a very cumbersome task if not impossible. 
Thus, in this paper, we derive a novel upper bound for the 
distribution of the sum of Double GG distributed RVs. Partic¬ 
ularly, a useful expression for the distribution of the product 
of Double GG distributed RVs is first obtained. Then, based 
on a well-known inequality between arithmetic and geometric 
means, a closed-form union upper bound for the distribution of 
the sum of Double GG distributed RVs is derived. This derived 
distribution is used to study the error rate performance of FSO 
systems with MISO/SIMO with EGC receivers employing 
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intensity modulation/direct detection (IM/DD) with on-off 
keying (OOK) over independent and not necessarily identically 
distributed (i.n.i.d.) Double GG turbulence channels. 

The rest of the paper is organized as follows: In Section II, 
the statistical characteristics of the Double GG distribution are 
provided, and an upper bound on the pdf of the sum of Double 
GG variates is presented. In Section III, we introduce the 
SIMO FSO system model. In Section IV, the BER expressions 
for SIMO FSO links are provided. In Section V, we present 
numerical results to confirm the accuracy of the derived 
expressions and demonstrate the advantages of employing 
spatial diversity over SISO links. Finally, Section VI concludes 
the paper. 

II. Definition and Statistical Characteristics 
A. Double GG distribution 

Let {I n }n = i be N statistically independent but not neces¬ 
sarily identically distributed Double GG RVs whose probabil¬ 
ity density function (pdf) follows (|T| at the top of the page 

H3 Eq. (5)]. 

In ([]}. G [.] is the Meijers G-functiorfl [{24} Eq.(9.301)], 
A (j-x) is defined as A (j; x) = x/j ,..., (x + j - 1 )/j , p n 
and q n are positive integer numbers that satisfy p n /q n = 
h,n/l 2 ,n, and trii^n > 0.5, i = 1,2, is a distribution shaping 
parameter. The distribution parameters 7 i n and il ir , of the 
Double GG model can be identified using the following 
equations 


The cumulative distribution function (cdf) of Double GG 
distribution can be derived from (|T]i as in © at the top of the 
next page HD- 


B. An Upper-Bound for the Distribution of the Sum of Double 
GG Variates 

Let us define a new R.V R, as the product of N Double 
GG R.Vs I n , i.e., 

JV 

r =n j «- (5 > 

n— 1 


The Double GG distribution considers irradiance fluctua¬ 
tions as the product of small-scale and large-scale fluctua¬ 
tions which are both governed by Generalized Gamma (GG) 
distributions, i.e. I n = U x , n U ytn , where U x , n and U Vtn are 
statistically independent arising respectively from large-scale 
and small scale turbulent eddies. Thus, R can be expressed as 
the product of 2 N GG R.Vs U n , i.e.. 


2N 
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n=1 


The pdf of R can be obtained using the statistical model 
proposed in (25ll as 


f R (r) = -A-Gj 


a? ,-1/3,0 


0,/3 


' j,OL 

— 

. W 

Ja (7l:2JV, mi :2 /v)_ 


(7) 


where £, w and J a (ji^N , tni-^N) are defined as 


Qi.n — 


r (m iin ) 


7i,i 


r (m i>n + l/7i,„) 

2 _ r (mi,„ + 2/7i >n )r 

<Tain “ T2 ( mi , n + 1/71,„ ) 

2 _ r (m 2 ,n + 2/72,™ )r (m 2 ,„) 
a V’ n - F2 (m 2 ,n + 1/72,n ) 


,71? i — 1? 2 

(2) 

*=( 


"n 

-1 

(3a) 

UJ = 

py faQ t 

a /71 

) 

-1 

(3b) 

l 

7~ \mai 



r (mi) 


Ja (7l:2AT, Wl:2Ar) — 


where cr^ n and a'y n are normalized variances of small and 
large scale irradiance fluctuations, respectively. The Double 
GG distribution accurately describes irradiance fluctuations 
over atmospheric channels under a wide range of turbulence 
conditions (weak to strong) H 2 - Furthermore, it is very 
generic and includes most commonly used fading models 
proposed in the literature as special cases such as Gamma- 
Gamma ( 7 = 1, = 1), Double-Weibull (rrii^ n = 1), 

and K channel ( 7 = 1 , 1 1^ n = 1 , TO 2 ,n = !)■ In addition, 
for the limiting case of 7 — >• 0 and to,.„ 00, Double GG 

pdf coincides with the log-normal pdf. 
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'Meijers G-function is a standard built-in function in mathematical 
software packages such as MATLAB, MAPLE and MATHEMATICA. If 
required, this function can be also expressed in terms of the generalized 
hypergeometric functions using (24] Eqs.(9.303-304)]. 
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is a positive integer with ki being also a positive integer. 
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IV. 

BER Performance 


The cdf of R can be derived from 0 as 
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geometric means, i.e. Am > Gm, with 
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a lower-bound for R.V Z dened as the sum of Double GG 
R.Vs, i.e.. 


N 
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can be obtained as 


Z > NR 1/n . (18) 

Considering Eqs. (IT4l) and (fl8l) . the cdf of Z is upper bounded 


as 
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By taking the first derivative of (fl9l > with respect to r, an upper 
bound for the pdf of Z can be obtained in closed-form as 

fz (r) < f* z (r) (20) 

where fz (r) is defined as 
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III. System Model 

We consider an FSO system employing IM/DD with OOK 
where the information signal is transmitted via one aperture 
and received by N apertures (i.e., SIMO) over the Double GG 
channel. We assume EGC receivers where the receiver adds 
the receiver branches. The received signal is then given by 
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= rjx ^2 In + v n , n = l,...,N 
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where x represents the information bits and can either be 0 or 
1, v n is the Additive White Gaussian noise (AWGN) at the 
n th receive aperture with zero mean and variance = (Vo/2 
, and p is the optical-to-electrical conversion coefficient. Here, 
I n is the normalized irradiance from the transmitter to the n lh 
receive aperture whose pdf follows 0. We should emphasize 
that the performance of SIMO under the assumption of equal 
gain combining is equivalent to that of MISO FSO links. 


A. Upper bound expression 

The optimum decision metric for OOK is given by ED 
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where r is the received signal vector. Following the same 
approach as m, ed, the conditional bit error probabilities 
are given by (see ED for details of derivation) 
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where 7 is the average electrical SNR obtained as 7 = if /Nq. 
Therefore, the average error rate can be expressed as 


PsiMO,ECG = 2 / h W erfc ( yi In 
j V n=1 


dl 
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where fi (I) is the joint pdf of vector I = (/i,/ 2 , ■ ■ • ,In)- 
The factor N is used to ensure that the sum of the N receive 
aperture areas is the same as the area of the receive aperture 
of the SISO link for a fair comparison. The integral expressed 
in d25l > does not yield a closed-form solution even for simpler 
turbulence distributions. However, an upper bound on (l25l) can 
be obtained by considering (l20l > and ( 12 1 1 ) as 
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The above integral can be evaluated in closed form by first 
expressing the erfc (.) in terms of the Meijer G-function 
presented in If26l eq. (11)] as 
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1 
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Then, a closed-form expression for (l26l > is obtained using 
Eq. (21)] as 
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The derived upper-bounded BER expression in (f28l> for SIMO 
FSO systems with ECG can be seen as a generalization of BER 
results over other atmospheric turbulence models. Specially, if 
we insert 7 * = 1 and H^ = 1 in (l28i l. we obtain an upper bound 
on BER expression over Gamma-Gamma channel. Setting 
mi = 1 in (|28]> . we obtain an upper bound for BER for Double 
Weibull channel. Similarly, for 7 , = 1, fi, = 1 and W 2 i = 1, 
an upper-bounded BER expression for K-channel is obtained. 


B. Diversity Gain and Asymptotic analysis 

Although Meijers G-function can be expressed in terms 
of more tractable generalized hypergeometric functions. (l28l > 
appears to be complex and the impact of the basic system and 
channel parameters on performance is not very clear. However, 
for large SNR values, the asymptotic behavior of the system 
performance is dominated by the behavior of the pdf near the 
origin, i.e. ( z ) at z —> 0 j27l . Thus, employing a series 
expansion corresponding to the Meijers G-function lE8l Eq. 
(07.34.06.0006.01)], (z) given in (EH can be approximated 

by a single polynomial term as 


/f(~) 


Na£, 


(uN aN ) 


aN\ mm 


f mi'y 1 

l a ’ 


n r fe 

i=i 

jAk 


Ck) i m 2JV72iv}-l 


(31) 


where Ck and c :] are defined as 

. JTO 171 TO 2 AT 72 An 

c fe = min - , f (32) 

t a a J 

c i e {A (“/ 72 ; TO2), • • •, 

A /a/ M\ • / TO l7l m-2Nl2N\ 

A ( /~f 2N ; m 2 N)} \ mm < - ,•••, -(33) 

la a ) 

It must be noted that (EH is only valid for independent and not 
identically distributed Double GG turbulence channels. Based 
on Eqs. (E 6 l > and (EH and at high SNRs, the average BER can 
be well approximated as 


-PsiMO.ECG 


r((l + Nac k )/2)t X 
2^c k (coN^) Ck y 
Uk 


c k ) 


(2N\ NaCk 


(34) 

Therefore, the diversity order of FSO links with N receive 
apertures employing equal gain combining is obtained as 
0.5 A min {771171 • • • , m 2 jv 72 jv}- 


V. Numerical Results 

In this section, we present analytical and Monte-Carlo 
simulation results using the previous mathematical analysis for 
the performance of SIMO FSO links employing EGC receivers 
over Double GG channels. The performance improvements 
over SISO systems are further quantified. Similar as in ll23l . 



Fig. 1. Average BER of EGC and SISO for plane wave assuming i.i.d. 
turbulent channel defined as channel b. 



Fig. 2. Average BER of EGC and SISO for spherical wave assuming i.i.d. 
turbulent channel defined as channel c. 


we consider the following four scenarios of atmospheric 
turbulence conditions reported in ifTTl 

« Channel a: Plane wave and moderate irradiance fluc¬ 
tuations with 71 = 2.1690, 72 = 0.8530, mi = 0.55, 

m 2 = 2.35, Oi = 1.5793, = 0.9671, p = 28 and 

q = 11 

• Channel b: Plane wave and strong irradiance fluctu¬ 
ations with 71 = 1.8621, 72 = 0.7638, m\ = 0.5, 

m 2 = 1.8, Hi = 1.5074, H 2 = 0.9280, p = 17 and 

q = 7. 

• Channel c: Spherical wave and moderate irradiance 
fluctuations with 71 = 0.9135, 72 = 1.4385, mi = 2.65, 
m 2 = 0.85, Hi = 0.9836 and H 2 = 1.1745, p = 7 and 























Fig. 3. Average BER of EGC over two i.n.i.d. atmospheric turbulence 
channels defined as channel a and channel b. 


q = 11. 

• Channel d: Spherical wave and strong irradiance fluc¬ 
tuations with 71 = 0.4205, 72 = 0.6643, mi = 3.2, 
m 2 = 2.8, fii = 0.8336 and Q 2 = 0.9224, p = 7 and 
Q = 11 - 

Figs. 1-2 demonstrate the average BER over i.i.d. channels 
defined as channel b and channel c, respectively. In order to 
verify the tightness of the bound, we present upper bound 
analytical results obtained through (l28l) along with the Monte- 
Carlo simulation of d25l) . As a benchmark, the average BER 
of SISO FSO link obtained through fl7l Eq. 24] is also 
included in these figures. As clearly seen from Figs. 1-2, the 
numerical results for the bounds are close to the equivalent 
simulated ones which represent the exact BER. For instance, 
at a target bit error rate of 1CF 5 , gaps between the exact and 
the upper bound curves are 5.2 dB and 6.6 dB respectively 
for N = 2 and 3 receive apertures employing EGC for 
channel b. Similarly, for channel c, at a BER of 10 , the 

gaps are respectively 3.8 dB and 4.3 dB for N = 2 and 3 
receive apertures employing EGC. This observation clearly 
demonstrates the accuracy of the proposed bound. It is also 
illustrated that the upper bound becomes tighter as the number 
of receive apertures decreases. This is expected as the upper 
bound and the exact BER curves coincide for N = 1. In 
addition, we observe that multiple receive apertures deploy¬ 
ment employing EGC signicantly improves the performance. 
Specially, at a target bit error rate of 10 -5 , we observe 
performance improvements of 46.8 dB and 66.8 dB for SIMO 
FSO links with N = 2 and 3 receive apertures with respect to 
the SISO transmission over channel b. Similarly, for channel 
c, at a BER of 10~ 5 , impressive performance improvements of 
51.1 dB and 63.9 dB are achieved for SIMO links with N = 2 
and 3 employing EGC compared to the SISO deployment. 



Fig. 4. Average BER of EGC over two i.n.i.d. atmospheric turbulence 
channels defined as channel c and channel d. 

Figs. 3-4 illustrate the BER performance of SIMO FSO 
links employing EGC receivers over non-identically dis¬ 
tributed (i.n.i.d.) Double GG channels. Similar to i.i.d. results, 
our upper bound closed-form expression yields a close match 
to simulation results. For example, at a BER of ICC 5 in SIMO 
links with N = 2 over i.n.i.d. channels a and b, the difference 
between the exact and the upper bound curves is 7 dB. For 
i.n.i.d. channels c and d and at a BER of 10 -5 , this gap is 5.8 
dB with N = 2. We also compare the performance of i.n.i.d. 
case with respect to i.i.d. case presented in Figs. 1-2. As an 
example, to achieve a BER of HP'’ in SIMO links over i.n.i.d. 
channels a and b, we need 8.5 dB less in comparison to i.i.d. 
case as channel a is less severe than channel b. Note that in Fig. 
1, we assume that both of the two channels between transmitter 
and receivers are described by channel b. On the other hand, 
to achieve a BER of 1CP 5 for SIMO links with N = 2 over 
i.n.i.d. channels c and d, we need 6.8 dB more in comparison 
to i.i.d channels as the channel d exhibits harsher conditions 
than channel c. Note that in Fig. 2, both of the channels 
between the transmitter and receivers are described by channel 
c. It can be further observed that asymptotic bounds on the 
BER become tighter at high enough SNR values confirming 
the accuracy and usefulness of the asymptotic expression given 
in 01. 

VI. Conclusions 

In this paper, we have derived a closed-form union upper 
bound for the pdf of the sum of Double GG distributed RVs. 
Using this bound, we have investigated the BER performance 
of FSO links with receive diversity employing equal gain 
combining over Double GG turbulence channels. An efficient 
and unified upper bound for the average BER of SIMO FSO 
systems with EGC receiver has been obtained which general¬ 
izes BER results over other atmospheric turbulence models as 

















special cases. Based on the asymptotical performance analysis, 
we have further derived diversity gains for SIMO FSO systems 
under consideration. We have presented BER performance 
based on analytical and numerical simulation results. Further 
comparisons between numerical and analytical results have 
confirmed the accuracy and usefulness of the derived results. 
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